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^'' ABSTRACT 

The gravitational lensing equations for convergence, potential, shear, and flexion are 
simple in polar coordinates and separate under a multipole expansion once the shear and 
flexion spinors are rotated into a "tangential" basis. We use this to investigate whether 
the useful monopole aperture-mass shear formulae generalize to all multipoles and to 
flexions. We re-derive the result of Schneider and Bartelmann that the shear multipole 
m at radius R is completely determined by the mass multipole at R, plus specific 
moments Q-™ and Qout OI " ^ ne m ass multipoles internal and external, respectively, to 
R. The m > multipoles are independent of Q ut- But in contrast to the monopole, the 
m < multipoles are independent of Qi n . These internal and external mass moments 



o 
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can be determined by shear (and/or flexion) data on the complementary portion of 
the plane, which has practical implications for lens modelling. We find that the ease 
of E/B separation in the monopole aperture moments does not generalize to m / 0: 
^ the internal monopole moment is the only non-local E/B discriminant available from 

lensing observations. We have also not found practical local E/B discriminants beyond 
the monopole, though they could exist. We show also that the use of weak-lensing data 
to constrain a constant shear term near a strong-lensing system is impractical without 
strong prior constraints on the neighboring mass distribution. 

Subject headings: gravitational lensing — methods: analytical 
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1. Introduction 

Weak gravitational lensing measurements of the shear 7 are often used to constrain the mass 
distributions in galaxies, cluster of galaxies, or even larger-scale objects. An exceptionally useful 
set of aperture mass formulae give non-parametric relations between the monopole moments of the 
lensing shear and the lensing mass. The aperture-mass formulae have these interesting aspects: 



1. A relation between the mean tange ntial shear co mponent 74 on a circle of radius R, and the 
mean convergence at and within R (|Kaiserlll995l ): 
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Recall that the convergence k is the surface mass density in units of the lensing critical density. 
The monopole of the tangent shear component at R is dependent upon the mass at or interior 
to R in a simple way, and is independent of the mass exterior to R. 



2. A rearrangement of the aperture-mass formula is (JFahlman et al.lll994l ) 



<«) 



<R 



drr (7 t ) r . 



(2) 



R 



This relation allows a model-independent determination of the mass monopole within R using 
on ly shear measures a t r > R. A generaliza tion to radially weighted aperture masses is given 
bv lKaiser et al.l (|1994l ) and ISchneiderl (|1996l ). which can allow the shear integral to have finite 
support. 

3. The shear monopole admits an instantaneous test for the presence of "B-mode" deflections. 
If we allow the lensing potential to be sourced by the normal scalar mass distribution ke plus 
a pseudo-scalar mass Kb, then we find that the tangent-shear formulae (pQ) applies to only 
the E-mode mass. The B-mode mass produces a monopole of the "skew" shear component 7 S 
rotated 45° from the tangent directional The skew shear is a direct test for B-mode sources 
within the aperture: 

(1s)r = {kb)<r - {hb)r- (3) 

Hence in real observations the monopole skew shear should be null. This is a special case of 
the general rul e that any E-mode mass measurement should be nulled when all shears are 



rotated by 45° (IStebbins et al. 



1996 



Luppino k, Kaiserlll997l ). 



In this paper we ask: can these three useful formulae relating monopole shear mome nts to 
monopole mass moments be extended to multipole moments? ISchneider &: Bartelmannl (119971 ) [SB97] 
offer a division of the shear multipoles into internal and external terms (their Appendix B), gener- 
alizing property (1). They further extend property (2), the ability to determine the mass moment 



The skew shear component is frequently designated by the misnomer "radial shear." 



from a closed-form integral of the shear, to the general multipole case. We offer here a simpler 
re-derivation of their results, extending them to the case where B-mode lensing may be present. 
In the process we also inquire whether property (3) can be extended: is there a simple test for 
-B-mode mass in the shear multipole signals? 

We fu rther ask: are the re equivalent properties for the higher-order lensing distortions, i.e. 



"flexions" ([Bacon et al.ll2006l )? The expected answer is yes, since these are sourced by the same 



two scalar degrees of freedom ke,b that produce the shear field. 

In $2] we derive a simple differential relation between the moments of convergence (mass), shear, 
and flexion, using a very compact notation for the standard lensing equations. This derivation will 
allow for both E- and B-mode source terms. In $3] we will examine the multipole generalizations of 
Equations ([I])-©. In SJH we give an application of these multipole formulae to a common problem 
in lens modelling: a galaxy-scale strong lens is embedded in a more extended group or cluster 
potential. We show how weak shear measurements could be used to constrain the cluster potential 
without assuming a particular geometry for the cluster mass. 



2. Lensing formulae in polar coordinates 

First we recast the familiar lensing equations into pola r coordinates. In t he flat-sky limit, 



lensing is compactly described by the differential operators of ICastro et al.l (J2005) 
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For an arbitrary deflection field (a x ,a y ) defined on the plane of the sky, we define a complex 
deflection a = a x + ia y . The deflection field can be decomposed into a curl-free E-mode part 
and divergence-free .B-mode part by defining a complex potential ip = ipE + iipB from the scalar 
potential ipE and pseudoscalar potential i/jb- 

a = dip. (6) 

The shear imposed on the background sources is given by the derivatives of the deflection: 

7 = 71 + il2 = -da. (7) 

Furthermore the complex convergence k = ke + i^B, which is the source term for the complex 
potential, can be written as 2k = da. Then the convergence and shear can be expressed in terms 
of the projected potential ip as 

2k = V 2 iP = dd^j (8) 

2 7 = ddi;. (9) 



From this we immediately derive the (|Kaiserlll995l ) relation between convergence and shear, which 
holds even when £>-modes are present: 

8k = &y. (10) 

The shear components defined with respect to the radius vector can be expressed as 

T = lt + i ls = - ie - 2ie . (11) 

Substituting this into the Kaiser relation (|1(J|) yields 
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where the subscripts after the comma denote differentiation, as usual. From Equation §§§, the 
shear can also be expressed as 



(13) 



Both of these equations clearly separate into radial and azimuthal parts under a multipole decom- 
position of the relevant quantities. For any complex quantity z we define the multipoles as 
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(1 + 6 m )ir 
The Kaiser relation can be rewritten for each multipole as 
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2.1. Flexions 



T he third derivatives of the lensing poten tial can be described by two complex-valued "flexion" 
fields (jCastro et al.1 120051 ; iBacon et al.1 12000 ): the spin-1 field T = ddd(f>/2, and the spin-3 field 
Q = ddd(j)/2. These satisfy 

F = dn, Bg = ddK. (18) 



In analogy with the shear, we define flexions in a tangential basis: 

F = e- ie T, G = e- 3ie g. 



(19) 



Casting Equation (|18p into polar coordinates, performing a multipole decomposition of F and G, 
then defining H = G — F, we obtain formulae for flexion multipoles in terms of convergence: 

or 
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Note the k dependence on the left-hand side is identical to Equation (|17p ; not surpringly, the flexion 
multipoles are very close to the shear multipoles. In particular 



(21) 



H (m) = 2 ( m + 2 ) T (m) 
r 

3. Multipole formulae 

3.1. Interior and Exterior Shears 

Several useful results may be obtained by integrating Equation (fT7|) by parts. First, we obtain 
a closed- form expression for the shear multipoles: 
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For m > we can assume Rf + V m )(i?i) -» and i?™ +2 r( m )(i?i) -► as i?i -» for any mass 
distribution which remains finite and differentiable at the origin. In this case 

rW(fi) = - /t M(fi) + 2( "' m + 2 1) f rdrr m ^ m) {r) (m > 0). (23) 

Since T^ m ' is completely determined by the mass distribution at and interior to R for m > 0, this 
is the desired generalization monopole formula dTJ. T(r) is, however, a complex quantity, so it is 
not fully specified by the m > multipoles. For m < a bounded mass distribution will have 
R™ +2 K ( m \R 2 ) -» and i^ +2 r( m )(i? 2 ) -»• as ^2 -> 00. We then obtain 

rM(fl) = - K ( m \R) - 2 ^^ f°°rdrr m K {m) (r) (m < 0). (24) 

« -Tt 

The negative-m multipoles of V are hence dependent only on mass at or exterior to R. This 
formula and the previous one fully specify the shear field. So the shear in a region R\ < r < R2 
is determined completely by the mass distribution in this region, plus these multipole moments of 
the mass interior and exterior to the annulus: 

Q-m\ R ) = I d 2 rr m e- ime K(r) (25) 

Jr<R 
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= TrR 2 (K) <R ((x-iy) m }, (27) 

Q^t(R) = [ d 2 rr- m e ime k(t) (28) 
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These definitions are normalized to agree with Eqns (B5) of SB97, who similarly demonstrate that 
shear at R depends upon interior and exterior masses only through these quantities. We have 
altered the phase conventions, however, in order to work successfully with the complex (E and B) 
convergence and potentials. The brackets in Equation (|2T[) indicate a mass- weighted average inside 
radius R. 

The present derivation shows that the division into interior and exterior shears holds even when 
there is an imaginary (S-mode) component to the potential and convergence. The shear components 
7i, and 7 S are always real- valued, as is k when there is no B-ioaode lensing: R^~ m > = K^ m >, etc. But 
this relation need not hold for T^ m > . It remains true, however, that the m > multipoles are 
produced by mass internal to R while m < are produced by external mass. 



Eqns (j23|) and (|25|) can be restated as 



rC")(r) = -^(r)+ (1 ^+^ +2 QL m) (r) (m > 0) (30) 
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For the flexions, Equation (|2ip makes it clear that H^ m > depends on local, internal, and external 
mass exactly as T^ m ' does. Equation (|20p shows that the F flexion depends only upon the local 
value of K^ m ' and its first derivative, independent of both the internal and external moments Q^ m '. 



3.2. Mass multipoles from Shear 

Multiplying Equation (fT7|) by r~ m , integrating by parts, and taking the upper integration limit 
to infinity yields 

^ m \R) = -T^ m \R) + 2{m + l)R m r drT (m \r) r - m - 2 (m > 0). (32) 
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Comparison to Equation (|23|) yields 

%n ' rdrr m K ( - m \r)=R 2m+2 r dr r - m - 2 T^ m \r) (m > 0). (33) 
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This is the desired generalization of the monopole aperture-mass formula ([2|). Multiplication by 
a weight function before the integration by parts would yield the full weighted aperture-multipole 
formulae of SB97. The mass interior to R affects the shear exterior to R only through the multipole 
moments Q-™ (R) ■ Here we see that these moments are completely recoverable from the shear field 
p("i) exterior to R. The Q^ are thus a complete description of the information that lensing data 
exterior to R can offer on the mass distribution interior to R. This holds even in the presence of 
B-mode lensing. 

Taking m < in this integration by parts yields an analogous formula by which the Q^(R) 
values may be determined from shear data at r < R: 

r\(~ m )(T>\ poo pR 

Vout W = / r drr m K ( - m \r) = R 2m+2 r dr r - m ~ 2 T {m \r) . (m < 0). (34) 
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When the shear data is available in a finite annulus Ri < r < R2, we can take differences of 
the above two formulae: 



/ d 2 r[ lt {r)+i ls {*)]r- m ~ 2 e-' ime = Rf m - 2 Q^\R 1 )-R^ 2m ' 2 Q^\R 2 ) (m ^) 
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3.3. Mass multipoles from Flexion 

Equation (|2ip implies that mass multipoles can be retrieved from the H^ m ' using the preceding 
shear formulae. Mass reconstruction from the F flexion differs: Equation (|20p integrates by parts 
to yield 

K M (R) = f " ^{r/R)- m F^(r) dr m > 

1 !o{r/R)~ m F( m \r)dr m< { ' 



3.4. E-B decomposition 

The real part of the monopole shear r'°) (R) at R determines the enclosed E-mode mass while 
the imaginary part gives the enclosed B-mode mass. This result does not generalize to other 
multipoles. Consider first the case where K,( m '(R) = 0, i.e. we are in a mass-free zone. Then 
for m > 0, r( m '(R) is fully specified by the complex number Q^'(R). But for m > 0, any 
chosen Q-™ amplitude and phase produced by an E-mode source ke can also be produced by a 
pseudo-mass source kb that is just the ke rotated about the origin by (90/m)°. There is hence 
no way to distinguish an internal m > E-mode mass distribution from an internal B-mode mass 
distribution. Similarly we can produce any desired Qout (-^0 with either ke or Kg source terms, so 



there is no test that can distinguish E-mode mass from B-mode mass distributions external to the 
shear measurement zone. These conclusions hold for flexion data as well as for shear data. 

The monopole turns out to be a special case: the E/B diagnosis is possible because the 
monopole E-mode and B-mode moments each have only 1 degree of freedom while the observable 
Q\J is complex. But for m 7^ 0, the E and B mass moments each have 2 degrees of freedom, so 
cannot be independently retrieved from a single Q. 

Thus if we have shear data on the R\ < r < R2 annulus, we have hope only of testing k for E 
and B components at m 7^ only within the annulus, not interior or exterior to it. Ideally this can 
be done by noting 

2 K P = K ^) + ^- m ) (38) 

2mf ] = K ^-R^ m \ (39) 

This can be combined with Equation (|20p . for example, to yield a pure-E quantity: 

2 (r/sgt + 4t + - 2 4 m) A) = r- m ^ (r^H) + r -|. (^j*-™)) . (40) 

Sending F — > iF gives a pure-B quantity. These equations are not practical null tests for B-modes, 
however, because they involve derivatives of F, which have divergent noise in the presence of shot 
noise from finite sampling. A practical null test for E/B modes in an annular region would require 
to an integral of F (or T or H ) over the annulus which could be approximated by a sum over source 
galaxies. We have not been able to derive such a form. 



4. Application to strong-lensing models 

In modeling a lensing system around a galaxy, one has strong-lensing constraints from multiply- 
imaged sources. The lens-mass model often contains a galaxy mass distribution K g (r), but it is 
essential in most cases to consider the influence of the larger-scale mass distribution on the system. 
Call this the "cluster" mass, which generates potential ip c { r )- On the assumption that the cluster 
mass has little structure on the scale of the strong-lensing system, V'c(r) can be approximated by 
a few terms of a Taylor expansion about the galax y center within some radius R\ that contains 



all of the strongly lensed features (JKochanekl Il99ll ) . The constant and linear terms of the Taylor 



expansion are immaterial to the strong-lensing model. The potential at r < R\, to cubic order in 
the Taylor expansion of the cluster, is 



V>(r) = (1 - k c ) 
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+ yr- 2 , (41) 



with V 2 ip g = 2n g . The strong- lens data produces a likelihood distribution over the (complex) 
parameters {7, o~, 5} and the galaxy-mass parameters. The mass-sheet degeneracy leaves k c uncon- 
strained by strong-lensing data. 



We now ask what additional constraints on these model parameters are available from the 
shear field at r > R\, We do not want to assume that the Taylor expansion offers an adequate 
description of the cluster mass at r > R\ , but from the previous discussion we know that only the 
multipole moments Q^(Ri) affect the parameters in the strong- lensing potential. In particular 
the terms of the form r rn e ± vm <t> [ n the Taylor expansion of (f) c can only be generated by mass outside 
R\ while the other terms can only be generated by mass inside R±. Specifically: 



The k c term is a monopole (constant) mass distribution, and Q 
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itR 2 k c . 



The 7 term is a constant shear, producing P 2 > = —27. From Equation (f3Tj) . (5„ utc (i?i) 

— irj. 



• The a term is a dipole mass distribution, k = Re(<rre ld ). This being the only dipole cluster 
mass kept in the expansion, we have Q inc = ivR\a/A. 

• The 5 term is an m = 3 external shear, r' -3 ** = — 26r, implying Q^ t C (R\) = —tt5/2. 

Eqns (|35p and (|36p can now be applied to a shear measurement that extends to radius R2 from 
the galaxy center. The multipole moments Q are split into galaxy and cluster contributions. Those 
from the galaxy are calculable from the parametric form adopted for K g . The cluster contributions 
at R\ are parameterized by the Taylor expansion coefficients as above. The cluster contributions 
at i?2 are formally unconstrained, but if R2 is large enough then these may be bounded by even a 
rough estimate of the total mass and extent of the cluster. We obtain: 
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+ (1 - K c ) -RlQ ou t,g( R l) + R 2Qout,g( R 2] 
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(42) 



(43) 



(44) 



(45) 



R lQout,q( R l) + R 2Qout,a( R 2 



Note that if the galaxy mass distribution has inversion symmetry about the coordinate origin, then 
all of the Q' 1 ' and Q( 3 ' moments of the galaxy vanish. In each equation, the left-hand side is an 
observable quantity and the right-hand side is a function of the galaxy parameters, the Taylor- 
expansion parameters, and some (presumably small) terms for the multipole moments to radius 
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i?2- Given a priori estimates of the R2 moments, each observation yields an additional constraint 
on the strong-lens model, with the mass-sheet degeneracy now broken. 

Were we to extend the cluster potential Taylor series to quartic terms, we would have 5 new 
degrees of freedom. The shear integrals for Q- m and Q ov [t wou ld constrain four of these. The 
ip oc r 4 term, however, would be an internal monopole mass distribution with r 2 radial dependence. 
It would contribute to Q in and would be degenerate with k c within the shear annulus. The 
degeneracy would have to be broken by the strong-lensing data. 

The shear data on the left-hand sides are affected by the intrinsic shape variation of the 
source galaxies, hence the constraints provided by each equation are not exact. As R2 — > 00, the 
negative powers of r in the integrands of Equations (I42p and (J43H lead to bounded shape-noise 
uncertainties on k c and a. Unfortunately the shape noise in the integrands of Equations (I44p and 
(|45p is divergent as the annulus grows outwards. Hence these expressions cannot be expected to 
provide useful constraints on the cluster quadrupole and octupole mass moments 7 and 5. 



4.1. (Constant) external shear 

Common practice in analysis of galaxy-scale strong-lensing systems is to limit modelling of 
external mass to a constant shear across the strong-lensing system, i.e. the 7 term of Equation (14ip . 
If we simplify Equation (|44|) by setting k c = and ignoring the shear induced by the galaxy mass, 
we find that the observable quantity on the left-hand side is just the mean shear inside the annulus, 
and 

7=j|^(^)-^^(7>ann. (46) 

One approach is to use a priori knowledge of the mass fluctuation spectrum to find a radius R2 

(2) 
beyond which we can expect the Q ^ t term to become negligibly small, < 0.01. Note that in this case 

the desired shear 7 is opposite to the mean shear in the annulus. Two problems arise however: first, 

it is not clear that any such radius exists, since the large-scale "cosmic shear" is typically ~ 0.01 

even before amplification by the R\jR\ factor in this term. Second, the shape-noise variance of 

the measured (i?| — R\){l)a,nn/R\ contribution will grow as (R2/R1) 2 for fixed R\, rendering the 

measurement uninteresting. It thus appears problematic to use weak-lensing information to infer 

the "external shear" in galaxy-scale lenses. 

In a different limit, one might assume that the mass within the R\ < r < R2 annulus has 
negligible quadrupole moment, perhaps because one does not see any galaxy groups or clusters 

(2) (2) 

projected within this annulus. In this case, Q Q ^ t (R2) = Q ^ t (R\) = —^7, and our estimate of the 
shear parameter becomes simply equal to the mean shear in the annulus 

7 = (7)ann- (47) 

In this case the shape noise on 7 decreases as R2 is increased and can become usefully small. 
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It thus appears that the use of weak-lensing information to constrain the external shear on 
galaxy lenses is practical only if one has a priori knowledge that an annulus in the vicinity of the 
lens is free of mass that would generate a significant shear on the system. We note that these 
problems are exacerbated for the higher-order external moments because of the positive powers of 
radius that appear within the shear integrals, e.g. Equation (|4"5|) . 



5. Conclusions 

Polar-coordinate expressions for the relations between convergence, shear, and flexion are 
seperable under multipole expansions once we rotate the shear and flexion spinors into "tangential" 
bases. Two well-known monopole aperture-mass properties are extensible to all m > multipoles: 
first, the shear multipole T^ m '(R) is determined solely by the convergence (mass) multipole kS 171 ' 
at or interior to radius R. The effect of the interior mass is fully described by the moments 
Q™ (R) oc J d 2 rr m K^ m >. Second, we find that the value of the interior mass moment Q^ (R) 
can be exactly recovered by an integral of the shear multipole F^ m ' from R to oo. 

The multipoles m < 0, however, behave oppositely to the monopole: the shear at R is deter- 
mined exclusively by the mass at or exterior to R. And the relevant exterior mass moment QoTt {R) 
can be determined by an integral of the shear interior to R. 

The tangential flexion component H = e~ Q — e J- behaves exactly as the tangetial shear 
r. In fact they differ only by a factor 2(m + 2)/r. The vector flexion component F depends purely 
on the local behavor of k, as is well known. 

The simple E/B decomposition of the monopole mass distribution does not generalize to 
m ^ 0. We show that shear or flexion data in an annulus R\ < r < Ri cannot discriminate 
between E-mode and B-mode sources outside this region — except for the monopole case. Shear or 
flexion data may be able to distinguish E from B sources inside the annulus, but we have not been 
able to derive a practical estimator which does this. 

The multipole formulae presented by SB97 and extended here find application in using weak- 
lensing data to constrain the large-scale characteristics of mass distributions in the vicinity of strong- 
lensing systems. Unfortunately a complete characterization of the exterior mass distributions is not 
practical because some of the aperture multipole formulae have divergent shape-noise behavior. In 
particular the estimation of the constant "external shear" term often found in strong-lens models is 
problematic without strong prior constraints on the neighboring mass distributions. We can expect 
the aperture-multipole formulae to find further use in generating model-independent measures of 
the shapes of dark-matter halos. 

This work is supported by grants AST-0607667 from the National Science Foundation, De- 
partment of Energy grant DOE-DE-FG02-95ER40893 and NASA BEFS-04-0014-0018. 
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